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ABSTRACT 

We study the time evolution of non-axisymmctric linear perturbations of a rotating 
magnetised neutron star, whose magnetic field is purely poloidal. The background stel- 
jy^ | lar configurations are generated self-consistently, with multipolar field configurations 

and allowing for distortions to the density distribution from rotational and magnetic 
O forces. The perturbations split into two symmetry classes, with perturbations in one 

class being dominated by an instability generic to poloidal fields, which is localised 
around the 'neutral line' where the background field vanishes. Rotation acts to re- 
duce the effect of this instability. Perturbations in the other symmetry class do not 
suffer this instability and in this case we are able to resolve Alfven oscillations, whose 
restoring force is the magnetic field; this is the first study of non-axisymmetric Alfven 
[ modes of a star with a poloidal field. We find no evidence that these modes form a 

J> ■ continuum. In a rotating magnetised star we find that there are no pure Alfven modes 

or pure inertial modes, but hybrids of these. We discuss the nature of magnetic insta- 
bilities and oscillations in magnetars and pulsars, finding the dominant Alfven mode 
from our simulations has a frequency comparable with observed magnetar QPOs. 

. Key words: 



1 INTRODUCTION 

Neutron stars are notable for the extreme strength of their magnetic fields, with surface fields reaching ~ 10 G for magnetars 
and interior fields perhaps being an order of magnitude stronger still. We expect many aspects of neutron star (NS) physics 
to be influenced by their magnetic fields, but we still have limited understanding of the actual structure of these fields; there 
are still open questions concerning their strength in the stellar interior, the relative proportions of poloidal and toroidal 
components and the possible effect of superconductivity (among others). One particular motivation for improving modelling 
of NS magnetic fields is the observation of quasi-periodic oscillations (QPOs) in the aftermath of giant flares from magnetars; 
these provide the first direct evidence of NS oscillations and give us a potential probe of the interior physics of these stars. 

One way to build up an improved understanding of NS magnetic fields is to explore the equilibria and dynamics of a 
simplifi ed NS model. We choos e to study a non-relativistic fluid star, but allowing for the effects of rotation and a magnetic 
field. In lLander fc Jones! we explored NS structures withm this model, generating self-consistent equilibrium solutions 



for stars with purely poloidal, purely toroidal and mixed poloidal-toroidal magnetic fields. To understand the stability and 
os cillations of these stars we perform t ime evolutions of perturbations, using our NS equilibria as background configurations. 



Lander. Jones fc Passamontil ( 2010h we studied the oscillation spectrum of NSs with purely toroidal magnetic fields, whilst 



In | , ., ,, v c c , 

Lander fc Jones! (|2010i ) presented results for toroidal-field instabilities. This paper discusses oscillations and instabilities of 



poloidal-field NSs, whilst in future we hope to complete this study of magnetised NSs by exploring mixed-field configurations. 

We begin by discussing the governing equations for our NS model, both for the background equilibria and for the 
time-evolution of the perturbations. We find that perturbations in one symmetry class are unstable in a region where the 
background field vanishes (consistent with the analytic work of, e.g., Wright (1973) and Markev fc Tavler 1 1973h ). with this 



unstable behaviour being reduced by the effect of rotation. We quantify the effect of magnetic field strength and rotation 
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on the instability's gro wth rate. Perturbations in the other symmetry class, by contrast, appear to evolve stably; this con- 
firms a prediction from Markev fc Tavler ( 19731 ) about the nature of the dominant poloidal- field instability. For these stable 
perturbations we are able to resolve many Alfven oscillation periods and extract their mode frequencies; this is the first 
study of non -axisymmetric Alfven modes in a poloidal-fie l d star, complementing a number of recent studi es on axisymmetric 
oscillations ( Sotani et al. 20081 : Sotani fc Kokkotas 20091 : Cerda-Duran et al. 20091 : Colaiuda et al. 2009h . We next turn to 
the oscillations of rotating magnetised stars, showing that these are a hybrid of inertial and Alfven modes. We conclude by 
discussing our work in the context of magnetars and pulsars and discuss the possible proportions of poloidal and toroidal 
components required for stability. 



2 GOVERNING EQUATIONS AND NUMERICS 

We begin by describing the equations governing our NS model, both for the axisymmetric stationary background configuration 
and the non-axi symmetric perturbations. A f u ller a ccount of the numerical methods used in the code and its performance 
may be found in Lander. Jones fc Passamontil ( 2010h : these are only described briefly here. 



We model a neutron star as a self-gravitating, rotating, magnetised polytropic fluid with perfect conductivity, in Newtonian 
gravity. We wish to study linear perturbations of this star; for this our governing equations consist of a set of stationary 
background equations and a set of equations describing the time evolution of the perturbations. The background configuration 
has a purely poloidal magnetic field Bo and may be (rigidly) rotating: 

= -VP - poV$ - PoO x (SI x r) + — (V x B„) x B , (1) 

V 2 $0 = 47TG/90, (2) 

Po = kp<, (3) 

where P is stellar pressure, p density, $ gravitational potential, G gravitational constant and SI angular velocity; 0-subscripts 
denote background quantities. Finally, we will take 7 = 2 throughout this study, as a rough approximation to a neutron star 
equation of state. 

Many studies of poloidal-field oscillations assume a dipolar field configuration, given by some simple analytic expression. In 
constrast to these, we solve for the field and fluid together, using a non-linear iterative procedure. The result is a self-consistent 
field configuration composed of a sum of different multipolar contributions. These higher multipoles are more significant in 
more distorted stars, providing small corrections to nonrotating and highly magnetised stars and becoming comparable with 
the dipolar field in rapidly rotating stars. Using the magnetic vector potential A defined through B = V x A, one may show 
that the magnetic vector potential (in spherical polar coordinates) satisfies the equation 

V 2 (A</, sin 4>) — — Kpor sin 6 sin <f>, (4) 

where k is a constant governing the strength of the field; a derivation of this is given in lLander fc Jones! \20o4 ). From the 
vector potential we obtain the background magnetic field configuration: 

B = V x (A^); (5) 

since B is poloidal it is described by a toroidal vector potential. These background equ ations are solved nume rically us- 
ing an iterative procedure t o find stationary equilibrium configurations, as detailed in lLander fc Jones! J2009) (see also 



Tomimura fc Eriguchil (2005)). Our method allows for distortions of the star due to both rotational and magnetic forces 



for the purely poloidal fields considered here, both of these forces will act to make the star oblate. 

For the perturbation equations, we work in the rotating frame of the background and write our equations in terms of 
the perturbed density Sp, the mass flux f = pov and a magnetic variable /3 = po<5B. We additionally make the Cowling 
approximation — neglecting the perturbed gravitational force — to avoid the computational expense of solving the perturbed 
Poisson equation. Our perturbations are then governed by seven equations: 

po^r = -jP V6p-2Slxi + ( (2 ~ 7)rP ° Vp - ^-(V x Bo) x B ) 5p+^-(\7xB )xf3+^-(\7xf3)xB -— — (Vp„x/3)xB , 

Ot \ po 4-7T / 47T 4-7T 47T/90 

(6) 

f-V.f, (7) 

§ = Vx (f xB )-^x (f xB ). (8) 
ot po 

The axisymmetry of the background allows us to decompose the perturbation equations in the azimuthal index m, thus 
reducing the three-dimensional system of equations to a 2D one. It is also beneficial to isolate perturbations of a particular 
m when discussing oscillation modes and instabilities. 
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2.1 Boundary conditions 

Instead of working with the spherical-polar radial coordinate r, we employ a coordinate x fitted to isopycnic surfaces, which 
for nonspherical stars is a function of r and 9. Doing this gives us a very simple set of boundary conditions at the stellar 
surface: 

t(x = R) = p(x = R) = , 8P(x = R) = 0. (9) 

Since we deal with ra^O perturbations, we should also enforce a zero-displacement condition at the centre: 

f(x = 0) = 0(x = 0) = , 8P(x = 0) = 0. (10) 

Perturbation variables are also zero at the pole, except for m = 1 where fe,f</>,/3g and are nonzero; in the case of these 
quantities, it is their ^-derivatives which vanish. 

Finally, it may be shown that perturbations of stars with purely poloidal background fields form two sets based on shared 
symmetry properties: Si = {f r ,f^,5p,/3g} and 5*2 = {fe, Pr, P<t>}- All the perturbations within a set have the same equatorial 
symmetry (say, even/symmetric), whilst all the perturbations in the other set have the opposite symmetry (which will then 
be odd/antisymmetric). We may enforce these symmetries as an additional set of boundary conditions, reducing our (2D) 
numerical domain to one quadrant of a disc. There are, therefore, two symmetry classes of solution: one class V + where 

elements of Si are even (with those in S2 odd) and one class V~ whe re .Si-elements are odd (with ^-ele me nts even). 

For more details on the boundary conditions summarised here, see Lander. Jones fe Passamonti ( 2010l ) or Lander fe Jonesl 



(2010). Finally, we note that whilst the poloidal background field vanishes on a 'neutral line' within the star (see section [3]), 



we do not impose vanishing of the perturbations there. 



2.2 Initial data 



Depending on the initial d ata chosen and the symmetry class specified (see above), either 'axial-led' or 'polar-led' perturbations 



(using the terminology of lLockitch fe Friedman! (| 19991 )) will be excited. In particular, axial-led perturbations are those in the 
class V~ and polar- led perturbations are in the V + class. For evolutions of polar- led perturbations we begin with an initial 
spherical-harmonic perturbation in the density: 8p ~ Yi m {9, 4>), whilst if we wish to excite axial-led perturbations we prescribe 
a 'magnetic' spherical harmonic perturbation in the velocity: v ~ VYi m x e r . 

Evolving perturbations on a background star with a poloidal field, we find that axial initial data strongly excites an 
unstable mode. Although some oscillatory behaviour can be seen too, the unstable mode dominates the evolutions and 
prevents us from extracting information about axial-led modes. With polar initial data, however, there is no evidence of an 
instability, and in this case we are able to resolve enough Alfven oscillations to study polar-led mode frequencies. For this 
reason, all results in the section on instabilities come from evolutions using axial initial data (and the corresponding symmetry 
class V~), whilst all results in the oscillations section use polar initial data (and the other symmetry class V + ). 



2.3 Numerics 

As described above, we need only evolve perturbations on one quadrant of a disc, by exploiting symmetries of the original 
3D system of equations. We employ a McCormack predictor-corrector algorithm to evolve the perturbation equations in 
time. We add an artificial viscosity term to our equations to damp out numerical instabilities that can emerge when using 
finite-difference routines, but we ensure that this is as small as possible — otherwise it could unnecessarily damp the physical 
hydromagnetic instabilities which we aim to study. We additionally employ a divergence-cleaning routine, si nce numerical 



error in the evolution of the perturbed magnetic field can result in an unphysical monopolar term: V • 5B 7^ (jDedner et al 
2002). We do not use artificial resistivity in this work. 



2.4 Dimensionless and physical quantities 

In our code we nondimensionalise all quantities by dividing them by the requisite combination of gravitational constant G, 
equatorial radius R and maximum stellar density p m . Some results presented in this paper are in terms of dimensionless 
quantities; this is either because they are frequently used in the literature or because there is no particular benefit to quoting 
them in physical units. Dimensionless quantities are denoted with a hat; for example Q — Q/^/Gp m . Note that dimensionless 
frequencies (i.e. the rotation Q, or mode frequency a) are quoted in terms of radians rather than cycles per unit code-time. 

In all cases where we have redimensionalised our results, this is done to the same 'canonical' neutron star, with a mass of 
1.4 solar masses and whose radius would be 10 km if it were non-rotating and unmagnetised (i.e. in hydrostatic equilibrium and 
spherical). For this canonical star the maximum density is 2.2 x 10 15 g cm~ 3 . Rotation frequencies may be redimensionalised 
using the approximate relation 

fi [rads -1 ] « 121000, (11) 
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Figure 1. Schematic field-line geometry for a poloidal magnetic field. Around the star's symmetry axis the field lines are 'open', in the 
sense that they extend outside the star; near the equatorial surface there is a region of field lines which close within the star. At the 
centre of these closed field lines is the neutral line (represented by the crosses in 2D, with the dashed line showing its circular form in 
the 3D star); along this line the magnetic field strength drops to zero. 



with the same relation holding for mode frequencies a. The corresponding rotation frequency v in cycles per second may be 
obtained using 

v [Hz] = - [r ^ ds ' 1] « 1900f2 (12) 

Z7T 

— again, the same relation may be used to obtain mode frequencies in Hz. 



3 INSTABILITIES OF PURELY POLOIDAL FIELDS 



The first indic atio ns that purely polo i dal fi elds could generically suffer hydromagnetic instabilities appeared in the studies of 
Wrightl jl973l ) and iMarkev fc Tavlerl \ 19731 ). Using coordinates adapted to the magnetic-field geometry, these authors found 
unstable perturbations in the vicinity of the 'neutral line', where the magnetic field vanishes; see figure [T] Instability was 
predicted to occur after approximately one Alfven crossing time — the time taken for a magnetic perturbation to travel 
around the star. 

The existenc e of un stable perturbations around the neutral line of a poloidal field is analogous to the toroidal-field case, 
for which iTavlerl (| 19731 ) showed the existence of an m = 1 instability localised around the magnetic fiel d's symmetry axis 



(the r egion where a toroidal field vanishes). A more general result for poloidal-field stability was gi ven by 



Van Assche et al 



Tavlerl (j 19731 ) 



( 19821 ). who showed that any poloidal field with closed field lines is unstable around the neutral line. IMarkev 
suggest that the open field lines, in contrast, may help reduce the effect of the instability. 

In contrast to the toroidal-field case, there does not seem to be a definite conclusion from these analytic studies about 
which values of azimuthal index m correspond to the most unstable modes; this is a consequence of them utilising field-line 
adapted coordinates rather than global spherical polars. One might, then, expect instab ilitie s to be present for a variety of m 
— and this was ind eed found in the num erical studies of Geppert fc Rheinhardt ( 20061 ) and Braithwaitel ( 20061 ) . 



Analytic work ( Pitts fc Tavler 19851 ) suggests that r otation may help to stabilise magnetic fields; this co nclusion was 
backe d up by the results of numerical evolutions reported in Geppert fc Rheinhardt ( 20061 ) . In contrast, the work of Braithwaitel 



( 20061 ) inds that rotation plays no stabilising role on poloidal fields; however, since this study models a star as spherical (with 
no distorting centrifugal force) it may be less relevant for describing fast-rotating, highly distorted stars. 

To summarise, the signatures of the poloidal-field instability described in earlier work are: localised unstable growth 
around the neutral line; onset after one Alfven crossing time; (probable) existence for a variety of m; (perhaps) stabilised by 
rotation. With these in mind, we now turn to our results for the behaviour of perturbations of a star with a poloidal field. 

To look for a poloidal-field instability, we begin by approximating the Alfven crossing time ta using volume- averaged 
background quantities: 

2R 



TA 



2R 



47T <p> 



(13) 



<c A > V B 2 

where R is the stellar radius, ca the Alfven speed and angle brackets denote volume averages. In dimensionless form, we find 



Oscillations and instabilities in NSs 5 



8AI 




5M 




Figure 2. m = 2 instability of a poloidal magnetic field in a 
nonrotating star. We plot the magnetic energy 8M against time 
t, both in dimensionless form, for three different grid resolutions. 
We see that the onset time for the instability is independent of 
resolution (appearing at around the expected value of ta ~ 77), 
and its growth rate converges, consistent with it being a physical 
instability. The results are for a star with average field strength 
B = 3.0 x 10 16 G. 



Figure 3. The stabilising effect of rotation on purely poloidal 
magnetic fields (for m = 2 perturbations). The magnetic energy 
is plotted against time for three different rotation rates. We see 
that increasing the rotation rate decreases the growth rate of the 
instability; i.e. the gradient of 8M is reduced in the regime where 
the instability dominates. As for the previous plot, each configu- 
ration has a field strength of B = 3.0 X 10 16 G. 
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Figure 5. The poloidal field configuration of the background star. 
The field strength decreases to zero around the neutral line, which 
is located at a dimensionless radius r/R ~ 0.8 from the centre. 



Figure 4. The magnitude of the perturbed magnetic field after 
the onset of instability, plotted with a logarithmic scale. The most 
unstable perturbations arc visible around the neutral line, where 
the background field goes to zero. This plot is for the m = 2 
instability, but very similar results emerge for m = 1 and m = 4 
perturbations. 



that ta ~ 77 for a star with a (redimensionalised) field strength of B = 3.0 x 10 16 G. In a plot of perturbed magnetic energy 
5M = J(SB 2 /8-k) dV against time, we would therefore expect to see an instability manifest itself around t « 77; this is 
confirmed in figure [2] To confirm that the origin of this instability is physical rather than numerical, figure [2] also shows that 
the result is similar for three different grid resolutions: low (16 r-points x 15 ^-points), medium (32 x 30) and high (64 x 60). 
The three growth rates converge at second order, the intended accuracy of the time-evolution code. 

In figure [3] we compare the behaviour of 5M for stars with a fixed field strength (3.0 x 10 16 G) but differing rotation 
rates, finding that rotation does act to slow the instability's growth rate. We will return to quantify the change in growth rate 
later in the section. 

We next turn to the prediction that an instability should be localised around the magnetic field's neutral line. From a 
global evolution, it is not straightforward to separate the behaviour of the instability from stable perturbations. We attempt to 
isolate the unstable perturbation by dividing the value of |<5B| at each point by its value close to the start of the evolution, in 
this way removing the shape of initial stable perturbations. We plot the resulting (dimensionless) quantity with a logarithmic 
scale. After the onset of instability, all plots we obtain are similar to that shown in figure SI for azimuthal indices m — 1,2 
and 4. Note that because our analysis is linear, the perturbed magnetic energy can grow indefinitely without diminishing 
the energy of the background configuration; in reality, nonlinear effects would become important after the initial period of 
unstable growth. 
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Figure 6. Left: the instability growth rate (in seconds -1 ) plotted against field strength B for a series of nonrotating stars; we see that 
the dependence is linear for these field strengths. Right: the effect of rotation on the growth rate of the m = 2 instability in a star 
with B = 3.0 X 10 16 . Note that the highest rotation rate shown is only ~ 58% of Keplerian velocity; beyond this value the code is less 
(numerically) stable. 



Comparing the shape of the unstable perturbation with the background field configuration — shown in figure [5] — we 
see that the instability is clearly dominant around t he neutral point. Furthe rmore, it appears to remain localised in this 
closed-field line region, fitting with the sugg estion of iMarkev fc Tavlerl l| 19731 ) that open field lines are likely to provide a 



stabilising influence. 

Having studied qualitative features of poloidal-field instabilities, we conclude this section with some quantitative results. 
As a measure of instability we define a growth rate £ by 

1 . A (&M\ 



where Mo is the magnetic energy of the background star. This simply measures the degree of exponential growth of the 
unstable mode and is easy to identify once the instability has come to dominate the behaviour of the system. In figure [6] we 
use C, to investigate the dependence of the instability on the background field strength and rotation rate of the star. From 
the left-hand plot we see that the growth rate scales linearly with field strength, with higher-m instabilities appearing to 
grow faster. In the right-hand plot of figure [6] we plot the m — 2 instability growth rate as a function of rotation frequency, 
up to Q ~ 0.58f2.R- . Rotation is seen to reduce the growth of this instability, but with finite-duration evolutions we cannot 
conclusively say that it will be removed altogether — however, we suggest that it may be stabilised in a real star by a 
combination of rotation and some additional physics (e.g. viscosity). A similar picture of rotational stabilisation seems to 
emerge for stars with Q > 0.6£Ik and for m — 1 and m = 4 instabilities, although evolutions in these cases appear to be more 
prone to late-time (numerical) instability. 

As mentioned in sections 12.11 and 12.21 each evolution of perturbations we perform uses one of two possible symmetry 
classes, dictating the equatorial symmetry properties of each perturbation variable. Perturbations in the V + symmetry class 
evolve stably for many Alfven timescales and we are able to study the corresponding 'polar-led' modes (see the following 
section). The generic instability of poloidal fields — explored in this section — manifests it self only in the behaviou r of V~ 
perturbations. We may understand this through a more detailed discussion of the study of Markev fc Tavler ( 19731 ). These 



authors work in a system of coordinates fitted to closed field lines, the geometry of which is a torus enclosing the neutral line. 
The y first 'open up' t h is tor us into a cylinder to enable comparison with known instabilities in a cylindrical field geometry 



(see IMarkev fc Tavlerl (|1973l ). figure 2). In the simplest case there are two main instabilities: the varicose instability (also 
known as the sausage instability) and the kink instability. Within a star the most unstable perturbations will be non-radial 
— those that do not involve work against the star's gravitational potential. Since the varicose instability does involve radial 
motion, the kink instability is likely to be the dominant one in stars with poloidal fields. In the spherical polar coordinates of 
the star, the kink instability involves motion of fluid elements in the ^-direction across the equatorial plane. Along this plane 
vg — for the symmetry class V + but not for the class V~ — and so unstable behaviour is only seen in the latter symmetry 
class. We believe this is the first numerical confirmation of this aspect of poloidal-field instabilities. 
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Figure 7. Polar-led m = 2 Alfven modes (a-modes) for a nonrotating star. We see that the mode frequency a varies linearly with B as 
expected, vanishing for unmagnetised stars (where the lowest-frequency mode is the /-mode). At the highest field strengths we see some 
possible deviation from the linear relation a oc B. We are unable to produce a similar plot for axial-led a-modes, as these oscillations 
were swamped by the dominant unstable mode. 



4 OSCILLATION MODES 

As described above, perturbations in the V~ symmetry class are dominated by unstable behaviour, preventing us from 
resolving this class of Alfven mode (a-mode) frequencies. In the V + class, however, the perturbations behave in a stable 
oscillatory manner; in this section we study their mode frequencies. All results reported in this section are for m — 2 modes 
— both for brevity and because these are easiest to resolve with our code. To check that these are representative of other 
azimuthal indices we performed evolutions for m = 1 and m = 4 oscillations too; the results suggested that there are no 
qualitative differences, but that m = 1 modes are slightly lower-frequency and m — 4 slightly higher-frequency than the 
m — 2 case. 

Before looking at res ults for oscillation mod e s, we discuss some terminology we will use to describe modes — this is based 
on notation employed by Lockitch fc Friedman ( 19991 ) to describe inertial modes (i- modes). The eigenfunctions of i- modes 



are — with the exception of the r-mode — a sum of spherical harmonic contributions Yj m and so cannot be labelled with a 
single index / (although they do have a single m). In all cases, however, the sums are dominated by the Yi m contributions 
from I = m up to some value Z = Zq ; the contributions beyond Iq are all far smaller. There may be several modes of the same 
m and lo, so these are enumerated with an additional index k and denoted mifc- Since we focus on m = 2 modes, we will 
suppress the m-index. Finally, we recall from section [2~T1 that we will h ave two classes of oscil l ation , corresponding to the two 
symmetry classes: these are axial-led {P ) and polar- led {P + ) modes. Lockitch fc Friedman ( 19991 ) showed that for axial- led 



modes the lowest contributing F; m to their eigenfunction (i.e. I = m) is axial and for polar-led modes the lowest term is polar. 



4.1 Modes of a nonrotating magnetised star 

In the nonrotating case, we find three peaks at lower frequency than the /-mode in the frequency spectrum; we identify these 
as magnetically-restored modes. We label these modes ai,a^ and as, based on the amplitude of their peaks, from strongest to 
weakest. We expect a-mode frequencies to be proportional to the Alfven speed ca = B/^/4-Kp and hence scale roughly linearly 
with B; this is borne out in our results, shown in figure [7] 

We can gain some understanding about the eigenfunction structure of these a-modes by comparison with the behaviour of 
the code for inertial modes (in unmagnetised stars). In this case we find that the lowest-Zo modes have the highest-amplitude 
peaks in frequency space — these modes are excited more strongly because of the finite resolution of the numerical grid, 
combined with the low-/ initial data we use. This provides us with a useful rough diagnostic to identify a-modes: we suggest 
that the eigenfunction of the strong peak ai contains lower Yi m contributions than the ai or az modes. 

In a perfectly conducting medium, like the model NS considered in this paper, magnetically-restored oscillations can 
occur in a continuous band of frequencies rath er than being d i screte global modes. This result was established by analytic 



work for an incompressible medium (see, e.g., iGoossens et al.l (|1985l )) and mor e recent numerica l work has suggested that 



the axisymmetric oscillations of compressible stars may form a continuum too ( Sotani et al. 20081 : Cerda-Duran et al. 20091 : 
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Figure 8. Polar-led m = 2 hybrid magneto-inertial modes, for a neutron star with field strength B = 6.0 X 10 G; the plot is clearer 
for this highly magnetised star than at lower values of B. When Q = there are three pure Alfven modes: oi, ai and a.3. Rotation splits 
each into co- and counter-rotating modes, which are seen to become known inertial modes for high rotation rates (with the possible 
exception of the highest-frequency mode; see text), although the large number of avoided crossings makes it difficult to track each mode 
individually. In these dimcnsionlcss units Keplcrian velocity Qk rj 0.72. 



Colaiuda et al .120091) . It i s known, however, that dissipative effects like viscosity and resistivity can act to remove the continuum 



Ireland et al 



(fl992T l) 



(e-g 

To test for a mode continuum, we look at the oscillation frequencies of perturbed quantities at different points within 
the star. If our system has discrete global modes we expect all these local oscillation frequencies to be equal; if there is a 
continuous mode spectrum then oscillation frequencies will be position-dependent. From our evolutions we find the former: 
mode frequencies at different points within the star are equal, within the resolution dictated by the length of our evolutions 
(of the order 1%). Whils t this appear s to contradict recent s t udies on magnetar QPQs that have found evidence for a mode 
continuum (among them lSotani et all (|2008h . lColaiuda et all (|2009T l and lCerda-Duran et al] (|2009h ). it should be emphasised 
that these studies are not quite comparable. For one, we study non-axisymmetric oscillations, rather than the m = modes 
in these earlier studies; in addition we have only looked at polar-led modes, in contrast with the axial oscilla t ions o f the other 
papers. Our work is, however, in agreement with the study of polar Alfven modes by Sotani fc Kokkotas ( 20091 ') . who also 
found a discrete oscillation spectrum. 



4.2 Modes of a rotating magnetised star 

Next we consider the mode spectrum of a r otating star with a poloidal magne tic field, but let us first recall our re- 
sults on toroidal-field oscillation modes from Lander. Jones fc Passamontil 1 201(J ). In this earlier work we found hybrid 



magneto-inertial modes, whose cha racter was Alfven-like for slow rotation and inertial-like in more rapidly-rotating stars 
( Lander. Jones fc Passamontil 2010h . so that their character depended on the ratio of magnetic to kinetic energy M/T — 



we expect to see this happen in the poloidal-field case too. In the toroidal-field case, we found two polar-led a-modes in 
the nonrotating case. After rotational splitting, one half of each a-mode appeared to become a zero-frequency mode in the 
M/T — > limit, whilst the other half of each mode became an inertial mode: 3 ii and 3 i2. This led us to classify the a-modes 
as analogues of their respective inertial modes, viz. 3 ai and 3 d2- 

In some respects, we find that polar-led oscillations of a poloidal-field star appear to differ qualitatively from the results 
for toroidal fields. In the non-rotating toroidal-field case we found two clear a-mode peaks, whereas with a poloidal field we 
find three — one of which is stronger than the others. Adding rotation splits these three modes into six — i.e. three corotating 
and three counter-rotating magneto-inertial modes; see figure [8] Close to Keplerian velocity Qk these all appear to become 
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known polar-led inert i al mo des: the two lo — 3 modes and the four lo = 5 modes of a rotating unmagnetised star (see, e.g. 
Lockitch fc Friedmanl l| 19991 )). 



The large number of avoided crossings (at which the character of two modes changes) present in figure [8] make it difficult 
for us to make conclusive statements about the eigenfunctions of a\,a2 and 0,3, because unlike the toroidal case we cannot 
easily track these modes from the Alfven-dominated to inertial-dominated regimes. However, since each a-mode splits into 
a co- and count er-rotating branch we e xpect them to become a corresponding pair of co- and counter-rotating i-modes. We 
also know from Lockitch fc Friedman (1999) that 3 ii, 5 ii, 5 i2 corotate with the star and 3 i2, 5 «3, 5 «4 are counter-rotating. 
Combining this information with the speculation that the strong peak a± is a lower-Zo mode than 02 and 03, we suggest that 
di becomes the pair of modes { 3 ii, 3 «2} as M/T — > 0, whilst 02 — > { 5 i2, 5 i3} and 03 — ¥ { 5 ii, 5 i4}. 

Finally, we note the slightly anomalous behaviour of the highest-frequency hybrid mode shown in figure [S] which begins 
at fl = as one branch of the et3-mode. This mode does not approach its apparent inertial counterpart 14 as closely as the 
other hybrid modes in the high-fi regime. Since the other five hybrid modes do convincingly become inertial modes, however, 
we suggest that this discrepancy is due to an avoided crossing with the corotating branch of the /-mode (not shown in the 
figure) at tl ^ 0.6. 



5 POLOIDAL-FIELD EFFECTS IN MAGNETARS AND PULSARS 



In the previous two sections we presented results for instabilities and oscillation mode frequencies at a variety of field strengths 
and rotation rates. Here we will apply these to two specific cases: a canonical model magnetar witlQ v — Hz and B = 10 16 
G and a model pulsar-NS with v = 1 Hz (a common, if slightly low value) and B = 10 14 G (assuming a fairly high surface 
field strength of 10 13 G); in both cases we assume that the volume-averaged magnetic field B is an order of magnitude greater 
than the surface field strength. 

We first consider the stability of poloidal fields in each model star. Although the magnetar field strength is 100 times 
that of the pulsar, a purely poloidal field would be unstable in both — however the instability growth rate, linear in B, 
would be 100 times slower in the pulsar, even if it were nonrotating. The additional effect of its rotation will further slow the 
growth of the pulsar's poloidal-field instability. For the magnetar we see from the left-hand plot of figure [6] that its growth 
rate £ « 300 s~ ; equivalently its e-folding timescale is around 3 ms. If the pulsar were nonrotating its e-folding timescale 
would be 0.3 s; rotation will increase this somewhat, though we cannot quantify this using our results. This is because our 
rotating background stars are constructed by specifying the oblateness of the star (which depends on the grid spacing) rather 
than its rotation frequency, so our slowest-rotating stellar models still rotate far faster than the 1 Hz of the model pulsar we 
wish to discuss here. 

Although a purely poloidal fie l d is u nstable for a number of azimuthal indices, the instability is highly localised; this 



was predicted bv lMarkev fc Tavlerl (| 19731 ) and our results (figures [4] and [5]) are in good agreement with this. As suggested by 
Wrightl (| 19731 ). a toroidal field which is fairly strong in the neighbourhood of the neutral line could remove any poloidal-field 
instabilities. Interesting, the toroidal field would not need to be globally strong to remove such instabilities. Based on this, 
we suggest that in a stable mixed-field configuration the maximum values of the poloidal and toroidal components may well 
be comparable but their respective energies need not be (the toroidal-field energy is likely to be considerably smaller, since 
the toroidal component need only occup y a small volume of the star). This suggests that the 'twisted - torus' magnetic-field 



configurations discussed by (for example) Yoshida et al. ( 20061 ). Lander fc Jones ( 20091 ) and Ciolfi et al. ( 20091 ) may be stable 
equilibrium solutions, despite having toroidal components whose energy is < 10% of the total magnetic energy. 

Since purely poloidal fields are generically unstable they are not candidates for long-lived magnetic field configurations 
in magnetars. For this reason, we cannot be confident that our results for polar-led Alfven oscilla tions should closely re s emble 
those of a real magnetar. In addition, it is unclear how many of the observed magnetar QPOs fsee lWatts fc Strohmaverl (|2007l ) 
for details of these) originate as Alfven modes of the interior, as opposed to elastic modes of the crust (for example). However, 
if future observations of magnetar oscillations do include frequencies best fitted by our results, it could indicate that these 
QPOs represent a-modes of dominantly poloidal-field configurations. 

The frequencies of the three polar-led a-modes ai, 02 and 03 exist in the ratio 7:4:11 (with a maximum discrepancy of 2%). 
For our model magnetar we would expect a strong Alfven QPO (corresponding to a\) at 83 Hz, with less significant QPOs 
at 47 Hz (02) and 130 Hz (as). It is interesting to note that our predicted value of 83 Hz is rather close to the two dominant 
magn etar QPO frequencies observed to date: 84 Hz for SGR 1900+14 and 92 Hz for SGR 1806-20 (|Watts fc Strohmaver 
20071 ). 

For our model pulsar we know that the oscillations will be hybrid magneto-inertial modes, but predicting whether they 
will be more Alfven-like or inertial-like is difficult — as before, this is because of finite resolution dictating the allowed rotation 



1 A typical rotation frequency for magnetars is v = 0.1 Hz (equivalently Q S3 5 X 10 5 in dimcnsionlcss units), but we will assume that 
such a small value is negligible when B ~ 10 16 G — as suggested by figures \E\ (right-hand plot) and [8] 
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rates for our background stars. To attempt to gauge the relative influences of il and B on the model pulsar's oscillation modes, 
let us return to figure [8] (for a 6 x 10 16 G star). In this figure the Alfven mode ai is seen to split into two hybrid modes, 
with frequencies ~ 10% different from their inertial counterpar t s at £ 1 = 0.31 [y = 600 Hz). We also know that hybrid-mode 
frequencies scale with M/T (see Lander. Jones fc Passamonti ( 2010l ) for more on this), and that M/T ~ B 2 /Q, 2 . For our 
model pulsar B is a factor of 600 smaller than in the plot and fl is also a factor of 600 smaller, so M/T is roughly equal to 
its value for the figure [8] star at Cl = 0.31; from this we suggest that the pulsar's oscillation modes will be dominantly inertial, 
with a magnetic correction of ~ 10%. 



6 DISCUSSION 



In this paper we have studied the behaviour of perturbations of a neutron star with a purely poloidal magnetic field. The 
background equilibrium configurations are generated self-consistently; the density distribution may be distorted by a combi- 
nation of magnetic and rotational forces and the magnetic field is multipolar and poloidal, generalising earlier work for dipole 
fields with a simple analytic form. 

A number of previous studies have shown that poloidal fields generically suffer instabilities localised around the 'neutral 
line' (where the background magnetic field vanishes); we confirm that these are also present in our global evolution of 
perturbations and for our particular poloidal field. The instability is present for a number of non-axisymmetric perturbations 
— we have studied it for azimuthal indices m = 1,2 and 4 — but only in one of the two symmetry classes of perturbation. 
We argue that this confirms a prediction from Markev fc Tavler (1973) about which modes are most unstable. Our results 



suggest that rotation has a stabilising effect on purely poloidal magnetic fields. 

For the stable symmetry class we are able to evolve perturbations for many Alfven oscillation periods and hence extract 
mode frequencies. These non-axisymmetric oscillations app ear to be global modes, in contrast with the continuous spectrum o f 
axisymmetric modes discussed recently by various authors ( Sotani et alj|200sl : Cerda-Duran et al. 20091 ; Colaiuda et al. 20091 ) . 
We find three Alfven modes, whose frequencies scale linearly with the field strength. These are rotationally split into pairs of 
modes, one of which co-rotates with the star and one which is counter-rotating. As the stellar rotation rate is increased (or 
the field strength is reduced) these become inertial modes; so in the case of a rotating magnetised star pure inertial modes 
and pure Alfven modes are replaced by hybrid magneto-inertial modes. 

There are a few qualitative diff erences between t he results presented her e and our earlier work on toroidal-field oscillations 
( Lander. Jones fc Passamontil 2010l ) and instabilities ( Lander fc Joneslboioh . Poloidal- field instabilities are not predominantly 



confined to m = 1 as toroidal-field instabilities are, but instead exist for a variety of m. The poloidal-instability growth rates 
we have studied (m = 1,2,4) are all lower than the m = 1 toroidal-instability growth rate, though of the same order of 
magnitude. It appear s from our result s that in the poloidal-field case the instability is more localised than the toroidal-field 
instability studied in Lander fc Jones ( 2010l ). Perhaps the most notable difference is that the poloidal-field instability only 
exists in one symmetry class of perturbations rather than both. In terms of oscillations, the character of hybrid magneto- 
inertial modes appears different if the background field is poloidal rather than toroidal. In the former case we find that each 
hybrid mode becomes inertial as the rotation rate is increased; in the later case some of the hybrid modes appeared to become 
zero-frequency modes in this limit. 

Our conclusion that purely poloidal fields are generically unstable is hardly surprising; however, it was not guaranteed 
that these localised instabilities would show up in evolutions of a discretised set of perturbations on a relatively coarse grid. 
Having shown that our code (at fairly low resolution) is adequate to study these unstable perturbations, we now discuss some 
other ti me evolutions o f magn et ar oscillations which have not found such instabilities and speculate on reasons for this. The 
work of Colaiuda et al. ( 20091 ) uses a coordinate system adapted to magnetic field lines in order to reduce the perturbation 
equations to a one-dimensional system. In these coordinates there is no coupling between different closed field lines, which 
together with a bo undary condit i on of not evolving perturbations on the neutral line allows for stable evolutions without 
numerical viscosity. Sotani et al. 1 2008h did not employ these field-line coordinates and had to use nu merical viscosity to 



suppr ess instabilities connected with evolving this intrinsically ID system in 2D. The related work of ICerda-Duran et al 



( 2009h does not mention instabilities and their plots of local QPO amplitudes within the star show that there is little or no 
growth immediately around the neutral line. Their work uses the anelastic approximation to suppress fluid modes, but it is 
not clear whether that should also remove poloidal-field instabilities. 

All of these studies are based on long-term stable evolutions of perturbations on a magnetised background star; further- 
more, they evolve axial perturbations, which our work indicates are particularly unstable for poloidal fields. The important 
difference between their work and ours appears to be that all three papers we have discussed specialise to axisymmetric 
(m = 0) perturbations, whereas we are currently only able to evolve perturbations with m ^ 1. Axisymmetric perturbations 
may simply be less unstable; evolving these perturbations with (for example) numerical viscosity may be sufficient to remove 
instabilities altogether. 

We have discussed some of the roles a poloidal field may play in both magnetars and pulsars. The unstable nature of 
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purely poloidal fields means that they are unlikely to be good models of any NS magnetic field, although the instability growth 
rate is considerably slower for pulsars than magnetars. We find that the dominant polar-led Alfven mode of a poloidal-field 
magnetar has a frequency of around 83 Hz, whilst pulsar oscillation modes will be hybrid magneto- inertial in nature. Although 
dominantly inertial, these pulsar modes may have magnetic corrections of ~ 10% for fairly typical field strengths and rotation 
rates. 

It has long been th o ught that a magnetic field may be stable with a suitable combination of toroidal and poloidal 
components (e.g. Tavler ( 1980h ). Since the poloidal-field instability only occurs close to the neutral line, we suggest that 
the addition of a toroidal component in this region alone may provide stabilisation. The toroidal component would have to 
be locally comparable in strength with the poloidal component, but need only occupy a small volume — in this case its 
contribution to the total magnetic energy could be quite modest. F r om this we have an indication that some recent models 
of mixed-field equilibria (| Yoshida et all bood ; [Lander fe Jone3l2009l : ICiolfi et aljlioogh may be stable, despite having only a 
small percentage of magnetic energy contained in the toroidal component. We hope to explore the stability and oscillation 
spectra of these mixed-field configurations in a future study. 



ACKNOWLEDGMENTS 

SKL acknowledges funding from a Mathematics Research Fellowship from Southampton University. This work was supported 
by STFC through grant number PP/E001025/1 and by CompStar, a Research Networking Programme of the European 
Science Foundation. We thank Kostas Kokkotas, Kostas Glampedakis and the anonymous referee for their helpful comments 
on this work. 



REFERENCES 

Braithwaite J., 2007, A&A, 469, 275 

Cerda-Duran P., Stergioulas N., Font J. A., 2009, MNRAS, 397, 1607 
Colaiuda A., Beyer H., Kokkotas K.D., 2009, MNRAS, 396, 1441 
Ciolfi R., Ferrari V., Gualtieri L., Pons J.A., 2009, MNRAS, 397, 913 

Dedner A., Kemm F., Kroner D., Munz C.-D., Schnitzer T., Wesenberg M., 2002, J. Comp. Phys., 175, 645 

Geppert U., Rheinhardt M., 2006, A&A, 456, 639 

Goossens M., Poedts S., Hermans D., 1985, Solar Phys., 102, 51 

Ireland R.C., Van Der Linden R.A.M., Hood A.W., Goossens M., 1992, Solar Phys., 142, 265 

Lander S.K., Jones D.L, 2009, MNRAS, 395, 2162 

Lander S.K., Jones D.L, Passamonti A., 2010, MNRAS, 405, 318 

Lander S.K., Jones D.L, 2010. larXiv: 1009.24531 

Lockitch K.H., Friedman J.L., 1999, ApJ, 521, 764 

Markey P., Tayler R.J., 1973, MNRAS, 163, 77 

Pitts E., Tayler R.J., 1985, MNRAS, 216, 139 

Sotani H., Kokkotas K.D., Stergioulas N., 2008, MNRAS, 385, L5 

Sotani H., Kokkotas K.D., 2009, MNRAS, 395, 1163 

Tayler R.J., 1973, MNRAS, 161, 365 

Tayler R.J., 1980, MNRAS, 191, 151 

Tomimura Y., Eriguchi Y., 2005, MNRAS, 359, 1117 

Van Assche W., Tayler R.J., Goossens M., 1982, A&A, 109, 166 

Watts A.L., Strohmayer T.E., 2007, Adv. Space Res., 40, 1446 

Yoshida S., Yoshida S., Eriguchi Y., 2006, ApJ, 651, 462 

Wright G.A.E., 1973, MNRAS, 162, 339 



